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Abstract 

Three-dimensional magnetic ordering transitions are studied theoretically in 
strongly anisotropic quantum magnets. An external magnetic field can drive 
quasi-one-dimensional subsystems with a spin gap into a gapless regime, thus 
inducing long-range three-dimensional magnetic ordering due to weak residual 
magnetic coupling between the subsystems. Compounds with higher spin 
degrees of freedom, such as N-leg spin-1/2 ladders, are shown to have cascades 
of ordering transitions. At high magnetic fields, zero-point fluctuations within 
the quasi-lD subsystems are suppressed, causing quantum corrections to the 
ordering temperature to be reduced. 
Compounds with strongly anisotropic crystal structures typically exhibit low-dimensional 
behavior at high temperatures. At low temperatures, the specific nature of their quantum 
fluctuations determines whether three-dimensional (3D) ordering jl] or other instabilities || 
occur, or whether there is no transition at all 0]. In particular, weakly coupled Heisenberg 
spin-1/2 chains are known to have a 3D magnetic ordering transition [|TJ or a spin-Peierls 
(SP) instability at low temperatures 0, if there is sufficiently strong coupling with low-lying 
phonon modes. On the other hand, compounds with an intrinsic spin gap, like e.g. weakly 
coupled integer-spin chains [f| or even-leg spin-1/2 Heisenberg ladders in a spin-liquid state, 
retain their one-dimensionality down to zero temperature || . The pure RVB nature of their 
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groundstate renders them inert to weak residual magnetic couplings between the quasi- ID 
subsystems. Thus 3D magnetic ordering and SP transitions are suppressed. 

In this second class of materials, an applied magnetic field, h, can decrease the singlet- 
triplet excitation gap of the quasi-lD subsystem, and eventually drive it into a gapless regime 
if the field exceeds a critical strength, h c \. A transition to a low-temperature ordered phase 
due to residual magnetic couplings becomes again possible in this partially polarized regime 
f||7|]. In this paper, we propose that 3D ordering as a result of the deconfinement of pairs 
of bound spinons by an external magnetic field can actually be realized in a wide variety of 
quasi- ID physical systems, including anisotropic spin chains, ladders, and SP compounds. 
Furthermore, a unified phenomenology for the magnetic phase diagram of these materials is 
presented, starting from an analysis of weakly coupled antiferromagnetic Heisenberg spin- 1/2 
chains (AFHC) with an easy-axis anisotropy. Interestingly, we observe that for compounds 
such as N-leg spin ladders with plateaus in their magnetization curves, m(h), || a cascade 
of N/2 ordering transitions for N even and (N+l)/2 transitions for N odd occurs at high 
magnetic fields. 

We are especially interested in the application of this magnetic field induced ordering 
transition to the material Cu 2 {C " 2 Hi 2 A ^ 2 ) 2 Cli (CuHpCl) ||. Typically, the spin gap in most 
of the ladder compounds known to date is too large to be overcome by presently available 
magnetic fields. However, this particular material has only a small spin gap of ~ 10.5 Kelvin 
which makes the interesting gapless regime experimentally accessible. It has recently been 
pointed out that CuHpCl may better be modeled by an ensemble of weakly coupled dimers 



than as an antiferromagnetic 2- leg ladder ||10|| . Whatever the precise magnetic structure 



may turn out to be, a magnetic field induced ordering transition can occur in all quasi-lD 
spin systems with a singlet-triplet excitation gap, including Ising-like chains, spin-Peierls 
chains, and ensembles of spin dimers. Other possible candidate materials with (relatively 



small) spin gaps include KCuCl% CuGeOs @, a' — NaV 2 0^, [[L2|], and the homologous 
series of cuprates Sr n _iCu n+ i0 2n [[J. 

Let us first consider a crystal of weakly coupled anisotropic antiferromagnetic spin-1/2 



, i quasi -ID subsystem 



FIG. 1. Schematic illustration of a strongly anisotropic 3D crystal, containing quasi-lD sub- 
systems, e.g. antiferromagnetic spin ladders. The coupling J' between the ID elements is small 
compared to the exchange constant J within the subsystems. Therefore the compound has quasi-lD 
properties at high temperatures (T 3> J'), while at low temperatures a ID to 3D ordering transition 
can occur if the ID subsystems are (driven) gapless. 



Heisenberg chains, described by the Hamiltonian 



H 1D = £ [j(Sf Sf +1 + SfS! +1 + &S?S* +1 ) - hSt 



(1) 



where J > is the antiferromagnetic exchange constant within the chains, A is an easy-axis 
anisotropy, and h is an applied external magnetic field. The chains are weakly coupled by 
H' = J' J2<i,j> S« ■ Sj with < J' < J. At zero magnetic field, the excitation spectrum of 
H 1D is gapless in the XY regime (A < 1) and massive in the Ising regime (A > 1), with a 
Kosterlitz-Thouless transition at the Heisenberg point (A = 1) fl3| . In the Ising regime, a 
finite critical magnetic field, h g , has to be overcome to completely soften the lowest triplet 
mode at wavevector tt, and to drive the system gapless. 

From a numerical solution of the Bethe Ansatz integral equations of Eq. 1 ||14|| , we 
have obtained the magnetic field dependent spinon velocity, u(h), the Luttinger exponent, 
K(h), the magnetization, m(h), and the radius of compactification, R(h), for the effective 
low-energy c=l conformal field theory. Magnetization curves, m(h), for an XY-like and an 
Ising-like anisotropic Heisenberg chain are shown in Fig. 2. In the dispersion curves of the 
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FIG. 2. Magnetization curve, m(h), of an antiferromagnetic spin-1/2 XXZ chain, (a) XY 
regime with Ising anisotropy: A = 0.5, and (b) Ising regime with A = 1.5. The insets show the 
dispersion law, u>(k), for the corresponding continuum models. Here the magnetic field acts as a 
chemical potential, rising from the bottom of the band at h c i to the top at h C 2- In the Ising regime, 
a spin gap opens around h=0, leading to a plateau in m(h). 

equivalent spinless-fermion model (insets of Fig. 2), the external magnetic field corresponds 
to a chemical potential. Below h c \ = — J(l + A) the cosine band, u)(k), of the spinless 
fermions is empty, and u(h c i) = (du/dk) h — 0. At the special field h — 0, the quasi- 
long-range ordered spin density wave within the chains becomes commensurate with the 
lattice, and Umklapp processes open up a spin gap for A > 1, leading to a plateau in the 
magnetization curve (Fig. 2 (b)) for —h g <h<h g , with 



where A = cosh(7) ||15| . Beyond h C 2 = J(l + A), the chains are fully polarized in the 
direction of the applied magnetic field. The particle-hole symmetry about h = is reflected 
by the shape of m(h) and the 3D ordering temperature, T c (h) (Fig. 3), discussed below. 

The low-temperature behavior of the susceptibility in the gapless regime is determined 
by the dominant low- frequency spinon modes at momentum q z = n: 
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FIG. 3. 3D ordering temperature (solid line) as a function of the applied magnetic field in a 
cubic crystal of weakly coupled antiferromagnetic spin-1/2 chains, (a) XY regime with A = 0.5, 
(b) Heisenberg point (A = 1.0), and Ising regime with A = 1.5. The dashed lines indicate the onset 
of the fluctuation region below which the 3D magnetic correlation length becomes comparable to 
the inter-chain spacing. For this plot, we have chosen J' /J = 1/16. 
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V8i^' AKJ u(l- 1/4*0 
where B(x,y) is Euler's beta function, and F(A) is a prefactor which strongly depends on 
the Ising anisotropy [JTSJ . Here, it has been assumed that the chains are parallel to the z-axis 
of the crystal. Furthermore, we have neglected higher-order logarithmic corrections which 
arise in a more rigorous treatment of the backscattering processes. 



Within the mean field approximation of Ref. JT^|, the low-temperature divergence in 



X (g 2 ,c<j = 0;T) drives a 3D ordering transition due to the weak inter-chain couplings, J': 



X 3D (q,^ = 0;T) 
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l + J'f(c L ) X 1D (q z ,uj = 0;Ty 
where /(q) is the crystal form factor which we here set to /(q) = —1 for simplicity (simple 
cubic lattice). The transition temperature, T c , is determined from the locus of divergence 



of x 3I) (q, uj = 0; T). In Fig. 3, the magnetic field dependence of T c is shown in the critical 
regime (A = 0.5), at the marginal point (A = 1), and in the massive region (A = 1.5). As 
the easy-axis anisotropy - equivalent to a nearest-neighbor repulsion in the spinless fermion 
picture - is increased, Umklapp scattering processes lead to a suppression of T c {h) around 
h = 0. In the Ising regime these processes become relevant, and a spin gap opens up, causing 
T c to vanish for \h\ < h g . 

To scrutinize the quality of this approach, we have calculated the Gaussian fluctuations 
which decrease the value of the transition temperature, T c (h), obtained from the mean 



field calculation [|18|]. This suppression is found to be rather small ((9(1%) ) because the 
corresponding Ginzburg-Landau parameter is proportional to the ratio J' /J <C 1. We 
thus conclude that the mean field treatment gives a quantitatively adequate description of 
the ordering transition. Furthermore, the onset of the fluctuation region below which the 
3D magnetic correlation length becomes comparable to the inter-chain lattice spacing is 
indicated by the dashed line in Fig. 3. This precursor region is quite large (C(~ 20%) ), 
indicating a sizable staggered magnetization, m s , in the ordered regime. Clearly, zero-point 
fluctuations are exponentially suppressed by J'. 

Let us now examine magnetic field induced 3D ordering transitions in quasi-lD systems 
with additional internal degrees of freedom, such as weakly coupled Heisenberg spin-S chains 
with S > 1/2, N-leg S=l/2 ladders with N>2, and spin-Peierls compounds. 

The low-energy Hamiltonian for a SP chain in a magnetic field is given by 

H ' SP = E [ J + S * • S *+i - hS i] > ( 5 ) 

i 

where 5 is the effective lattice distortion due to the coupling between spinon and phonon 
degrees of freedom. There are three well-known regimes in this system: (i) for \h\ < h c i, it 
is in a spin-liquid phase with a spin gap A 5G ~ 5 2/3 \ log5|- 1/2 PJ for sufficiently small 5's, 
and h c i = O(Asg)', (ii) for h c \ < \h\ < h c2 the system is a gapless spin-density wave with 
a field-dependent modulation; (iii) for \h\ > h c2 it is fully polarized in the direction of the 
applied magnetic field. 
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In the gapless, partially polarized, region (ii), we can map the low-energy spectrum 
of this system onto an effective XY-like Heisenberg chain (Eq. f) with the parameters: 
J e ff = J(l — S)/2, A e ff = 1/2, and h e ff = h — J(5 + 35)/ 4. This mapping is valid for 5 ^ 1. 
The lower and upper critical fields of the gapless regime are h c \ = J(l + 35)/ 2 and h C 2 = 
2 J |]2U|| . Considering weak antiferromagnetic 3D couplings between the SP chains, H' = 
J'J2<i,j> Si • Sj, the qualitative magnetic field dependence of the 3D ordering temperature, 
T c (h), is thus found to be the same as in Fig. 3(a). 

As seen in the above example, weakly coupled anisotropic Heisenberg chains in a field 
can serve as an effective theory for other quasi-lD quantum magnets with richer internal 
structure. Let us now extend this idea to weakly coupled N-leg spin- 1/2 Heisenberg ladders 
with the Hamiltonian 

H N = J\\ ^2 S«,r ' Sj, r + J±^2 Si,T ' St,r', ^i,Ti (6) 

<-* I i,r 

where i and j enumerate the rungs, r, r' label the legs, and the sum marked by <-> (\) 
runs over nearest neighbors along legs (rungs). Ladders with even and odd number of legs 
are known to have quite distinct features at zero field (h = 0): at low energies, odd-leg 
systems can be mapped onto effective gapless antiferromagnetic Heisenberg chains with 
longer-range interactions In even-leg systems, relevant inter-band scattering processes 



open up a spin gap for any finite positive inter-chain coupling, J± > [22| This is the 
reason why approximations based on strong coupling anisotropies, i.e. expansions in J\\/J± 
Pljm, give a correct qualitative picture, extending even beyond the isotropic (J\\ = J±) 



regime P,p2|. Furthermore, many physical ladder materials show coupling anisotropies 
within the ladder complex, e.g. a recent structural analysis of the candidate vanadate ladder 
material NaVa 2 5 suggests a strong rung-coupling anisotropy of J\\/J± ~ 13/75 [pl| , and 
for the proposed 2-leg cuprate ladder, CuHpCl, the anisotropy ratio was estimated to be 
J\\/J± w 1/5 J25J. 

It has recently been shown that N-leg Heisenberg spin- 1/2 ladders can have up to N/2 
plateaus in their magnetization curve, m(h), if N is even, and (N+l)/2 plateaus if N is odd 



TABLE I. Parameters of the effective low-energy model for the gapless regions of N-leg spin- 1/2 
ladders in a magnetic field, h. The effective magnetic field is given by h e ff = h — h c (0)J± — ChJ\\- 
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(0.302,0.2743, 


(0.3926,1.1506, 


(0.3432,0.2888, 




1.1917,1) 


0.1962,0.1875) 


1.6577,1.9239) 


0.5555,0.8125) 



S. These plateaus are related to gaps between the spin multiplets in the excitation spectra. 
The partially polarized regions between the plateaus are gapless, and can be described by 
effective massless c=l conformal field theories, corresponding to a spin- 1/2 AFHC Hamilto- 
nian with parameters {Jeff-, A e //, heff)- In table 1, the numerical values of these parameter 
sets, obtained from a strong rung expansion (J\\/J± <C 1) are listed for N = 2, 8 legs. 

Applying the mean field approach described above to a crystal of weakly coupled N-leg 
ladders, we now observe a cascade of N/2 ((N+l)/2) 3D ordering transitions for quasi-lD 
ladder subsystems with an even (odd) number of legs, shown in Fig. 4. In the case of weakly 
coupled even-leg ladders, the first transition is driven by the formation of a quasi-long-range 
ordered SDW of triplets along the ladder direction, commensurate with the magnetic field, 
h > h c \. The following transition (for A^ > 2) is driven by a SDW of quintuplets, etc. . 
Depending on the ratio J\\/J±, these phases of different multiplet polarization may overlap, 
and mixed phases can occur (see e.g. Fig. 4 (d) around h = 1.5 J± and 1.8JjJ. The 
resulting 3D ordering temperature does not vanish completely in this case, but has minima 
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FIG. 4. 3D ordering transition temperatures of N-leg spin-1/2 Heisenberg ladders as a function 
of an external magnetic field. Cascades of transitions are observed for N > 2, driven by ID SDW's 
of spin-S multiplets on the ladders. For this plot, we have chosen an anisotropy ratio J\\/J± = 1/5 
and residual inter-ladder couplings J'/J\\ = 1/16. 

at particular magnetic fields where the number of the lower multiplet excitations equals 
the number of the next-higher multiplet excitations. By analogy, a cascade in T c (h) should 
occur in compounds containing weakly coupled integer-spin Heisenberg spin chains. 

Odd-leg spin-1/2 ladders - as well as half-odd- integer-spin Heisenberg chains - also have a 
sequence of ordering transitions, with the only difference that the onset of the first transition 
occurs already at h — (Figs. 4 (b) and (c)). The dip features in T c (h), due to strong 
Umklapp scattering (A e jf < 1) become less pronounced at larger magnetic fields, i.e. A e jj — > 
0. This is expected because the effective spin degrees of freedom of the higher multiplets are 
more "classical" as S — > oo, and quantum fluctuations are suppressed at high fields. Finally, 
as the width of the ladder subsystems is increased, the gaps between the various multiplets 
disappear, and a quasi-continuous T c (/i)-curve emerges - as it is expected for the limiting 
case (N = oo) of weakly coupled 2D Heisenberg planes. 

In summary, we have examined 3D magnetic ordering transitions in quasi-lD materials. 
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In the presence of a spin gap, an external magnetic field can drive the ID subsystems gapless, 
thus inducing long-range 3D magnetic ordering due to residual magnetic coupling between 
the subsystems. For compounds with higher spin degrees of freedom, such as N-leg spin- 1/2 
ladders, cascades of ordering transitions are predicted to occur, reflecting the mesoscopic 
character of these materials. At high magnetic fields, quantum fluctuations in the quasi- ID 
subsystems are suppressed, causing the "dip feature" in T c (h) - due to Umklapp processes - 
to disappear for the transitions at higher fields. 

We wish to thank G. Bickers, C. Broholm, E. Dagotto, A. Honecker, D. Scalapino, and 
A. Tsvelik for useful discussions, and acknowledge the Zumberge Foundation for financial 
support. S.W. is partially supported by DAAD grant No. HSP III,D/98/11174. 
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